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Abstract 

A new precise measurement of \V c b\ and of the branching ratio 
BR(5° — > D* + £^U() has been performed using a sample of about 5000 semilep- 
tonic decays B° -> D* + £~i7 e , selected by the DELPHI detector at LEP I by 
tagging the soft pion from D* + — > D°tt + . The results are: 

= (39.0 ± 1.5 (stat.) (syst. exp.) ± 1.3 (syst. th.)) x 10~ 3 
BR(B° -> D*+£-z7^) = (4.70 ± 0.13 (stat.) toll (syst. exp.))% 

The analytic dependences of the differential cross-section and of the Isgur Wise 
form factor as functions of the variable w = Vb° ■ vd* have also been obtained 
by unfolding the experimental resolution. 
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1 Introduction 

In the framework of the Standard Model, the mixing between quarks of different 
flavours is described by the Cabibbo-Kobayashi-Maskawa (CKM) matrix. Its elements 
are not predicted by the theory, apart from the constraints due to the requirement of 
unitarity. 

A precise measurement of \V c t,\, the element corresponding to the beauty to charm 
quark transitions, constrains the parameters which describe the process of CP violation 
for mesons^ 0J. As a result of progress in the phenomenological description of heavy 
flavour semileptonic decays, \V c b\ is determined with small theoretical uncertainty, from 
either the inclusive process b — > c£~ug, or from an analysis of the form factors in the 
decay B% — > D* + £~i7 £ . The present measurement is based on the second approach 0. 

The decay rate for the last process is proportional to \V c b\ 2 and to the hadron matrix 
elements describing the transition from a B® to a D* + meson. In the limit of very heavy 
quarks (m^c >> Aqcd ~ 200 MeV/c 2 ), the amplitude is proportional to a single form 
factor jF(u> ), where w is the scalar product of the B® and D* + four- velocities. It is equal 
to the D* + Lorentz 7 factor in the B® rest frame. When w = 1, the D* + is produced at 
rest in the B® rest frame: as a consequence of Heavy Flavour symmetry, the normalisation 
•F(l) = 1 is expected. Corrections to this prediction due to perturbative QCD have been 
computed up to second order JJ. The effect of finite b, c quark masses has been calculated 
in the framework of the Heavy Quark Effective Theory [Q. The value J-"(l) = 0.91 ± 0.03, 
as determined in ||, was used in this analysis. This result is consistent with the value 
•F(l) = 0.88 ± 0.05, derived in || on the basis of the study presented in ]7|], and with a 
more recent computation based on lattice QCD which yielded .F(l) = 0.93 ± 0.03 0. 

The measurement of the decay rate at w — 1 would therefore determine \V c b\ with 
small theoretical uncertainty. Due to phase space suppression, this quantity is deter- 
mined from the extrapolation to 1 of the differential decay rate dT/dw , where T{w ) is 
parametrised according to several different functional forms P,|9]JiI| (see also discussion 
below). Results based on this approach have been reported by the ARGUS [ I I] and 



CLEO [|i"2||i"3]] collaborations operating at the T(4S) resonance, and by ALEPH [ 14 , 13 ] 



and OPAL |00 at LER The present paper updates the previous DELPHI result of 



reference ||18|| . Identification of D* + mesons is based on the tagging of the soft pion (n*) 
from the decay D* + — > D°tt + , the method referred to as "inclusive analysis" in reference 
18fl . As compared to this previous work, the following improvements were obtained: 



the resolution on w was improved by a factor of about 1.5 by applying the algo- 
rithm of inclusive secondary vertex reconstruction developed for B® lifetime [19| and 



oscillation [[2(J measurements; 

the full available statistics was analysed, thereby increasing the sample by more than 
a factor two; 

the most recent parametrisation [|1^] of T{w ) was used to extrapolate the experi- 
mental data to w =1; 

a more precise determination of the BR(6 — > i~ v £ D* + X) was used to compute the 
fraction of events in the sample due to non-resonant D* + tt production, or to the 
intermediate production of higher excited charm states which then decay into a D* + ; 
all these states will be called D** in the following. 



1 Charge conjugated states are always implied; lepton (£) means either an electron or a muon, unless the contrary is 
explicitly stated. 
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2 The DELPHI detector 



The DELPHI detector has been described in detail elsewhere ||21|| . Charged particle 
tracking through the uniform axial magnetic field (B = 1.23 T), secondary vertex re- 
construction and lepton identification are important in this analysis: they will be briefly 
described in the following. 

The detector elements used for tracking are the Vertex Detector (VD), the Inner De- 
tector (ID), the Time Projection Chamber (TPC), the Outer Detector (OD) in the barrel 
and the Forward Chambers in the endcap regions. The average momentum resolution for 
high momentum (p) charged particles in the polar angle region between 30° and 150° is 
a{p)/p = 0.0006 p (GeV/c) 0. 



The VD, consisting of 3 cylindrical layers of silicon detectors (radii 6, 8 and 11 cm), 
provides up to 3 hits per track (or more in small overlapping regions) in the polar angle 
range 43° < < 137°. In the original design the VD provided only two-dimensional in- 
formation in the R<p plane, orthogonal to the beam direction. Since the 1994 data taking, 
an upgraded detector with full three-dimensional point reconstruction was installed. In 
the R(j) plane the spatial resolution of the VD is about 8 pm. per point. Tracks from 
charged particles are extrapolated back to the beam collision point with a resolution of 



20 2 + 65 2 /p 2 \_ /^ m ; where p± is the momentum of the particle in the R(j) plane. The 
resolution on the z coordinate depends on z and is on average slightly worse than that 
in R<f). The primary vertex of the e + e~ interaction was reconstructed on an event-by- 
event basis using a beam spot constraint. The position of the primary vertex could be 
determined in this way with an average precision of about 40 pm (slightly dependent 
on the flavour of the primary quark-antiquark pair) in the plane transverse to the beam 
direction. Secondary vertices from B semileptonic decays were reconstructed with high ef- 
ficiency employing the algorithm described in reference ||19|| . The decay length resolution 
for the present analysis was about 400 pm. 

Leptons were identified among all the charged particles of momentum 2 < p < 30 
GeV/c. To allow the reconstruction of the B® decay point only particles with at least 
one hit in the VD were considered as lepton candidates. 

Electron identification was based on a neural network algorithm, optimally combining 
the information from the ionisation signal in the TPC, from the energy release in the 
electromagnetic calorimeters, and, for tracks with momentum below 3 GeV/c, from the 
Ring Imaging CHerenkov counters (RICH). A level of tagging providing about 75% ef- 
ficiency within the calorimeter acceptance was chosen. The probability for a hadron to 
fake an electron was about 1%. Electrons from photon conversions are mainly produced 
in the outer ID wall and in the inner TPC frame. About 80% of them were removed with 
negligible loss of signal by reconstructing their materialisation vertex. 

Muons were selected by matching the track reconstructed in the tracking system to 
the track elements provided by the barrel and forward muon chambers. The efficiency 
was about 80% for about 1% probability of hadron mis-identification. 

The experimental efficiencies and hadron mis-identification probabilities were mea- 
sured year by year using dedicated samples of leptons and hadrons independently tagged 
and the simulation was tuned consequently. 



3 Hadronic Event Selection and Simulation 



Charged particles were required to have a momentum in the range 0.25 < p < 45 GeV/ c, 
a relative error on the momentum measurement less than 100%, a distance of closest ap- 
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proach to the interaction point less than 10 cm in Rcj) and 25 cm along z, and a polar 
angle such that |cos# |< 0.937. Electromagnetic showers not associated to tracks were 
required to be well contained within the calorimeter acceptance and to have an energy 
release greater than 0.5 (0.3) GeV in the barrel (forward) electromagnetic calorimeter. 
Only hadronic showers with an energy release greater than 1 GeV and not associated to 
tracks from charged particles were accepted as neutral hadrons. 

The following selection was applied to the detector operating conditions: the TPC was 
required to be fully efficient, and at least 95% of the electromagnetic calorimeters and 
90% of the muon chambers had to be active. Hadronic Z decays were selected with 95% 
efficiency and negligible background by using standard cuts (see reference [[2l|1 ). 

Each event was divided into two opposite hemispheres by a plane orthogonal to the 
thrust axis. To ensure that the event was well contained inside the fiducial volume of the 
detector the polar angle of the thrust axis of the event had to satisfy the requirement 
|cos6> | < 0.95. Charged and neutral particles were clustered into jets by using the LUCLUS 



22| algorithm with default resolution parameter djoin = 2.5 GeV/c. 



About three million events were selected from the full LEP I data sets. The JETSET 



7.3 Parton Shower |22| program was used to generate hadronic Z decays, which were 
followed through the detailed detector simulation DELSIM [23| and finally processed by 
the same analysis chain as the real data. A sample of about seven million Z — > qq events 
was used. To increase the statistical significance of the simulation, an additional sample 
of about 2.2 million Z — > bb was analysed, equivalent to about ten million hadronic Z 
decays. Details of the Z samples used are given in Table [I]. 



Year 


real data 


simulated 
Z -> qq 


simulated 
Z ^ bb 


1992+1993 
1994+1995 
Total 


1203982 
1832082 
3036064 


2012615 
5190586 
7203201 


922764 
1321384 
2244148 



Table 1: Available number of events. In 1992 and 1993 only two-dimensional vertex 
reconstruction was available. 



4 The D* + t v sample 

4.1 Event Selection 

Only events containing at least one lepton candidate were considered further. The 
transverse momentum of the lepton relative to the jet it belonged to, pf, was computed 
after removing the lepton from the jet. The cut p[ > 1 GeV/c was imposed to reduce 
the background. 

A charm hadron candidate was reconstructed from all the particles in the jet containing 
the lepton, except the lepton itself, by means of the iterative algorithm described in detail 
in reference [|l9j . Small clusters were first formed out of the charged particles and, when 



possible, a decay vertex was computed for each cluster. The charm candidate so obtained 
was intersected with the lepton trajectory to provide the B% secondary vertex. In the 
case where only one charged particle with hits in the VD belonged to the cluster, its 
intersection with the lepton track was computed. The cluster associated to the secondary 



4 



vertex with the largest statistical significance Sc (defined as the distance from the primary 
vertex divided by its error; in years 1992 and 1993 only the projected distance onto the 
R(f) plane was considered) was kept as a seed. All other charged and neutral particles in 
the jet were ordered by decreasing values of their pseudo-rapidity relative to the cluster 
direction, and added to it provided the mass of the system did not exceed 2.2 GeV/c 2 . 
The charm three-momentum was obtained from the sum of all the particles assigned to 
the cluster. The charm trajectory was evaluated again and was finally intersected with 
the lepton track to obtain the B® decay point. To improve background rejection and the 
resolution on w (see below), events with significance iS,c<4.5 were rejected. 

The 7r* candidate (the pion from D* + decay) was searched for among all particles in 
the jet with charge opposite to that of the lepton. If the candidate belonged to the charm 
cluster, the D° four-momentum was computed after removing the it* from the cluster 
and imposing the D° mass. To increase efficiency, particles classified as fragmentation 
products were also considered as tc* candidates. The D° was then identified with the 
charm cluster, constrained to the D° mass. D* + production was finally tagged based 
on the mass difference Am= M^o^* - M^o (see Figure [I]). All events with Am<0.165 
GeV/c 2 were used for the analysis. 

4.2 Event Kinematics 

The variable w (= % • vjj*+) can be expressed as: 

d 

w = (M| + M 2 D , + - q 2 )/(2M §0 M D * + ) 
where q 2 is obtained from the and D* + four-momenta as: 



Q 2 = (Pbo ~Pd*+) 2 

d 

The D* + energy, polar and azimuthal angles, and the energy of the B® meson were 
determined as in reference The resolution obtained in the simulation was: 



<r{E § o)/E §0 = 10% 

d d 

a(E D *+)/E D *+ = 12% 
(t(9d*+) = 18 mrad 
a((f)D*+) = 21 mrad 

The B® direction was evaluated using two estimators: 

• the direction obtained by inverting the vector sum of all the particles in the event 
except the ones assigned to the B® . This procedure, already used in reference [18 



exploits three-momentum conservation in the event. The resolution depends on the 
hermeticity of the detector, but can also be spoiled whenever another semileptonic 
decay takes place in the event; 

the direction of the vector joining the primary and the secondary vertex: the res- 
olution achieved depends on the distance between the two vertices, improving for 
higher values. This approach was not used in the inclusive analysis of reference 



Using the simulation, the resolution was parametrised on the basis of the missing energy 
in the hemisphere opposite to the B® for the first estimator, and as a function of the 
reconstructed decay distance of the B® for the second. The B® meson direction was 
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then obtained as the average of the two, weighted by the inverse of their error squared. 
When the difference between the two values was greater than three times its error, the 
direction nearer to the D* + l~ system was chosen. In the years 1992 and 1993 only the 
first estimator could be used to determine the B® polar angle. The resolution function 
obtained could be parametrised by a Breit-Wigner distribution, with half width at half 
maximum: 

r(0g- o )/2 = 12 mrad 

d 

r(0 B - o )/2 = 12 mrad (1994- 1995) 

d 

r(0 B -o)/2 = 24 mrad (1992 - 1993) 

d 

The resulting w resolution function is shown by the dots in Figure 0. The RMS width 
of the core of the distribution is approximately the same for all data sets (a(w ) = 0.125), 
but larger tails are present in the 1992-1993 sample due to the poorer 9 measurement. The 
RMS width corresponds to about 25% of the allowed kinematic range (1 < w < 1.504). 
Due to resolution effects (17.9 ± 0.4)% ((32.9±0.6)%) of the events of the 94-95 (92-93) 
data set lay outside that range. 

The squared recoil mass fi 2 was also determined on the basis of the event kinematics. 
It is defined as: 

/i 2 = M% + M 2 D , +i - - 2P B - • P D * n - , (1) 

d d 

where M^ ^ D * +l -y P§o^ D * +i -^ are the mass and four momenta of the B® meson and 

D* + t~ system respectively. In the decay process B\ — > D* + £~i7 e , jj? represents the 
square of the mass of the neutrino, and should be zero. In the case of background pro- 
cesses, due to the emission of additional particles other than £~ , v and D* + , it is usually 
greater than zero. The cut jj?< 2 GeV 2 /c 4 was applied to reduce the D** contamina- 
tion. The square recoil mass was also used to improve the w resolution: the constraint 
fi 2 = M 2 (= 0) was imposed on equation (|T]), which was then inverted to improve the 
determination of the B® polar angle 8^ - A second order equation was obtained: the 

d 

resulting ambiguity was solved by choosing the solution nearer to the previous deter- 
mination. When the resolving discriminant was negative, it was forced to zero. This 
procedure improved the precision on the determination of w both for 1992-1993 and for 
1994-1995 data samples, reducing the amount of B% — > D* + £~i7 e decays outside the al- 
lowed kinematic range to (4.8±0.1%). The shaded area in Figure ^| shows the w resolution 
finally obtained. 

4.3 Sample Composition 

A set of J\f% = 10232 events was finally selected. One contribution to the background 
was the combinatorial component, due to random association of a hadron and a lepton. 
Another was the resonant one, due to the association of the lepton to a true tc* produced 
by processes different from B® — > D* + £~v e . 

The combinatorial background was determined from the real data, by applying the 
previous selection to all candidates in the jet containing the lepton and having the same 
charge as the lepton (wrong sign sample). A few events in this sample were in fact due 
to resonant processes, when either the lepton from D° semileptonic decay or else a fake 
hadron with the same charge as a true n* was selected. Their respective amount was 
computed from the simulation as 63 ± 31 events and 76 ± 38 events. The total yield 
of 139 ± 49 events was subtracted from the wrong sign data set. This sample was then 
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normalised to the right sign sample by counting the events situated in the side band 
interval 0.225 GeV/c 2 < Am < 0.3 GeV/c 2 , where the fraction of events due to genuine 
7i* was negligible. The normalisation factor was 1.288 ± 0.012 (stat.)± 0.021 (syst.) and 
the corresponding number of combinatorial events in the mass interval selected for the 
signal was N CO mb. = 3737 ± 70 (stat.)±75(syst.). The systematic error consisted of three 
contributions. The first one (±54 events) was computed by applying the same procedure 
to the simulated data, after having removed all events containing a genuine 7r*, in order 
to verify that the particles with wrong charge correlation reproduce the actual combi- 
natorial background. The difference between the right and wrong charge samples, after 
normalisation, was 44 ± 54 events. The second contribution (±49) was due to the sub- 
traction of the small amount of resonant events in the wrong charge sample (see above). 
The residual contribution, due to leakage of n* events into the side band, was negligible. 

The total amount of D* + was then ftf D *+ = 6495 ± 123(stat.) ± 75(syst.). 

The following processes contributed to the resonant background: fake leptons ran- 
domly associated with a tt*, b decays to a D* + with another heavy flavour decaying 
semileptonically, b — > D* + X c /r~X (followed by X c /r~ — > £~Y), and production of D** 
in b semileptonic decays. The contribution from all these sources was determined from the 
simulation, using the most recent measurement of the relevant branching ratios, which 
are reported in Table (|. 

Hadrons faking a lepton can combine with a D* + produced either from bb or cc de- 
cays of the Z (the contribution from gluon splitting to D* + is negligible). Their total 
amount was computed by determining independently the probability for a hadron to fake 
a lepton, known with about ±5% relative precision, and the product of branching ratios 
BR(Z -> bb{cc)) x BR(6(c) — > D* 



The rate for the B® — > t~v t D* + decay was obtained from the measurement of the 
inclusive BR (6 — > t~i7 t X c ) =(2.6 ± 0.4%) [|25|| , multiplied by the probability that the 
charm state (X c ) hadronises to a D* + . This last number was estimated as (50±10)% 
from the fraction of B® semileptonic decays with a D* + in the final state ||. 

The fraction of inclusive double charm decays b — > D* + X c was determined from charm 
counting measurements as suggested in reference f26fl . The error on the signal included 
the uncertainty in the inclusive semileptonic decay X c — > l~V(X. 

The main contribution to the resonant background is due to the intermediate pro- 
duction of D** states. In the following it will be assumed that the D**^D* + X de- 
cay is saturated by single particle production (namely, D**° — > D* + n~, D** + — > D* + ~ J ' 



TV 



D** s — > D* + K°), a hypothesis consistent with the conclusions of reference [27]. Therefore 
the rate for D** background production at LEP can be expressed as: 

b** = BR(b -> £-u e D**) x BR(D** -> D* + X) 
= f u ■ BR{B~ r vtD* + <K-) ± 
/„ ■ BR(5° i-U.D^n ) + 
f 8 ■ BR(B S -> ru e D* + K°) 

where the parameters f u , fd, f s express the probability that a b quark hadronises into a 
B~ , B° and B s meson respectively (the production of D* + from semileptonic decays 
is neglected). Their values are reported in Table ^ as computed in reference |J. The 
relation f u = f d is also imposed. 

By assuming that the partial semileptonic widths are the same for all b hadrons, the 
following relations are also derived: 



BR(B° d -> £~i7 e D* + 7r ) = BR(B~ i~i7 e D* + 7r- 



1 r B» 
X - X 



2 T B - 
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BR(B S -> £-u e D* + K°) = BR(B~ -> £~i7 £ D* + Tr-) x(-a) x (2) 

4 t b - 

where the factor | in the first relation accounts for isospin invariance, the factor | in 
the second one is derived from SU(3) flavour symmetry. The parameter a = 0.75±0.25 
is introduced to account for a possible violation of the SU(3) symmetry. The branching 
ratio BR(B- -> (ri7 e D* + Ti-) was determined by the DELPHI J27| and ALEPH J28J col- 
laborations by looking for an additional charged pion coming from the B~ decay vertex 
in a sample of exclusively reconstructed £~ D* + events. The two measurements are in 
good agreement and provide the average value: 

BR(b -> B- -> £-u e D* + n-) = (4.76 ± 0.78) x 1(T 3 (LEP) (3) 



The ARGUS collaboration JTTJ] has determined the fraction of D** in their sample of 
B — > D* + £~X events from a fit to the fi 2 spectrum. Using the same model assumptions 
as in their paper, the value 

BR(B% -> £~i7 e D* + ir ) = (6.2 ± 1.9 (exp.) ± 0.6 (model)) x 10~ 3 (4) 

is derived, which corresponds to 

BR(b -»• B~ -> r^D* + 7r~) = (5.3 ± 1.7) x 10" 3 (ARGUS) (5) 

Equations (0) (||) and (|5|) are finally combined to provide: 

b** = (0.76 ±0.11 ±0.03 ±0.02)% (6) 

where the first error is experimental, the second is due to the error on the B hadron 
production fractions (mostly / s ) and the third comes from the variation of the parameter 
a in the range 0.5 - 1. 

The events generated in the simulation were rescaled to the branching ratios deter- 
mined previously and were then processed through the same analysis chain as the real 
data. This allowed the determination of the composition of the selected D* + sample 
which is reported in Table 0. 



Source 


Amount 


Data 


10232 


Combinatorial 


3737±70 


B° d -> D*+ti7 t 


54±3 


b -> D*+X c 


56±3 


fake leptons 


250±8 


b -> D**£~i7 e 


1469±10 


B° d -> D* + £-U e 


4666±130 



Table 2: Expected composition of the sample used for the analysis. Only the statistical 
errors are reported. 



5 Determination of \V c b 
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5.1 Parametrisation of the Decay Width 

The expected number of signal events can be expressed as a function of w by the 
relation: 

dM/dw = 4N Z R b f d BR(tt*) e(w) dT/dw, 

dT/dw = — M^ + (M B -o - M D , + f Vw^l {w + l) 2 (7) 



x I Vrh ? T 2 {w 



d 

4w 1 — 2wr + r 2 
1 + w (1 — r) 2 



The factor 4 accounts for the fact that a B® can be produced in either hemisphere, and 
that both electrons and muons were used; Nz is the number of hadronic events, Rb the 
fraction of hadronic Z decays to a bb pair, fd the probability for a b quark to hadronise 
into a B® meson, BR(7r*) the branching ratio for the decay D* + — > D°tt + , t^ the B^ 

d 

lifetime and r is the ratio of meson masses, r = Mr>*+ / M^ - The values employed for 

d 

these parameters, as determined by independent measurements, are reported in Table 
f|. The quantity e(w), the product of the acceptance and of the reconstruction efficiency 
(which exhibits a slight dependence on w ), was determined on the basis of the tuned 
simulation. 

The analytical expression of the form factor J-"(w ) is unknown. Because of the small 
w range allowed by phase space, earlier analyses used a Taylor series expansion limited 
to second order: 

F(w) = ^(1) (1 + p%{\ -w) + c(l - wf + 0(1 - wf) (8) 

Except for .F(l), theory does not predict the values of the coefficients, which must be 
determined experimentally. First measurements of \ V c b\ were performed assuming a linear 



expansion, i.e. neglecting second order terms []TT| , |T2| , P~4l, [T8|| . Basic considerations derived 
from the requirements of analyticity and positivity of the QCD functions describing the 
local currents predict however that a positive curvature coefficient c should be expected, 
which must be related to the "radius" of the heavy meson p?r(see reference 0) by the 
relation: 

c = 0.66,4-0.11 (9) 
Results exploiting this analyticity bound have been derived by the ALEPH and OPAL 



collaborations (see reference |]I5| , |I6|| ). 

An improved parametrisation was subsequently proposed (see reference ||). It ac- 
counts for higher order terms, so reducing to ±2% (according to the authors) the relative 
uncertainty on \V c b\ due to the form factor parametrisation. In this approach, the four- 
velocity product is first mapped onto the variable z, defined as: 

The form factors are then computed by continuing z in the complex plane, where it is 
bound to lie within the unit circle. The form factors are then expanded as a power series 
of z while analyticity bounds and dispersion relations are employed to express terms up 
to order three as functions of the first order coefficient. The resulting expression is rather 
complicated 0. 
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An equivalent approach was applied in reference [10 where the form factors are ex- 



pressed instead as a function of w . In this case a novel function A\(w ) was introduced, 
connected to ^(w ) by the following relation: 



F 2 {w) x 
Al(w) x 



1 + 



4w 1 — 2wr + r 2 
1 + w (1 — r) 2 



(10) 



1 — 2wr + r 
! (1-r) 



2 r 



1 



n 2' 



where Ri(w ) and R 2 (w ) are ratios of axial and vector form factors; their analytical 
expressions can be found in reference [lOR. The following parametrisation, depending 



only on a single unknown parameter , was obtained for A\\ 



w 



A 1 {w) = At{l) 1 - 8p 2 Ai z(w) + (53p 2 Ai - 15)z(w) 



(231pi 



91)z(w) 



(11) 



where the relation between z and w was given previously. It should be noted that in the 
limit w —>■ 1, Ai(w ) — > ), so that Al(1) ~ ^"(l)- Experimental data were fitted 
using this last parametrisation and results were also obtained with the other forms for 
the sake of comparison. 

Results using this new parametrisation have also been published by the OPAL col- 
laboration in [17], and presented by the CLEO collaboration at the XXX International 



Conference on High Energy Physics El 



5.2 Parametrisation of the D** Spectrum 

The D** sample is composed of four different charm orbital excitations, two narrow 
resonances (Di,D$), with a measured width of about 25 MeV/c 2 p5| , and two broad 



states Pij-Do)- According to HQET, their width should be about 200 MeV/c 2 ; the 
CLEO collaboration has reported preliminary evidence of the D\ state with a mass of 
2461± 50 MeV and width = 290±100 MeV j2jj. Non-resonant D* + n states may also 
be present and contribute to the sample: it will be assumed in the following that their 
behaviour is included in that of the broad states. 

The differential decay width for the decay processes b — > D**£ui has not been measured 
and must be taken from theory. Heavy Quark Effective Theory predicts that, in the limit 
of infinite b,c masses, the rate near zero recoil ( w — 1) should be suppressed by an extra 
factor (w 2 — 1) compared to the B% — ► D* + £~i7 e decay rate. Several computations of 



the relevant form factors have been performed in this approximation (see reference ||30 
and references therein). However such models predict a high production rate for the D% 
states, which is incompatible with present experimental information (see discussion in 
reference f5I|l ). 

The effects of C(AQC£)/m c ) and 0(a s ) corrections have been computed in reference 
|32|| (referred to as LLSW model in the following). Decay rates and differential decay 
widths are computed assuming two different expansions, ("A" and "B" schemes) and 
the results are compared to the prediction obtained in the infinite mass approximation 
("Aoo" and "Boo"). A few parameters are not predicted by that model, but are varied 
within a reasonable range. When including finite c mass corrections, the D\ production 
rate decreases and is consistent with present experimental limits, while the D** rate near 
zero recoil is increased. The form factors for the broad states are computed from those of 
the narrow states assuming a non-relativistic constituent quark model with a spin-orbit 
independent potential. The model predicts that the global production rate for the broad 
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states should be about equal to that of the D\. Experimental results show that narrow 
resonances account for about (35 ± 15)% of b — > D**£u£ decays ||, in fair agreement with 
that prediction. It should be noted that in the "Boo" scheme the rate of broad states is 
1.65 times larger than the D\ one. However, the prediction for the D\ state is too high. 

A calculation based on a relativistic quark model ( "EFG" model) reduces the number 
of unknown parameters of the model (see reference ||3"5|]). In such a case, however, broad 
states account for only about 25% of b — > D**£p e decays. 

The following prescription, elaborated by the LEP \V c b\ working group, was applied 
to determine \V c t,\. Among all the possible expansions of the LLSW model, only those 
consistent with the experimental constraint from the ratio of the D\ to D\ production 
rates were considered. This removes the Aoo and Boo models. Each of the remaining 
models was then used in turn; input parameters were varied one at a time, while keeping 
all the others fixed at their central value. The allowed range for each parameter was 
once again determined from the D\ over D\ rate. The average of the two extreme | V^,| 
values so obtained was used as the measurement result, while half their difference was 
considered as the systematic error. 



5.3 Fit to the Experimental Data 

Real and simulated data were collected in several w bins. A minimum \ 2 fit was then 
performed comparing the numbers of observed and expected events in each bin. The nor- 
malisation of the background was determined as explained previously. The shape of the 
w distribution for the combinatorial background was obtained from the wrong charge real 
data events. Simulated D** spectra were corrected as described in the previous section 
and the spectra for all the other background sources were taken from the simulation. The 
contribution from the signal was obtained at each step of the minimisation by properly 
weighting each generated event surviving the selection; for a given value of w the weight 
was equal to the ratio between the value taken by the fitting function and the one of 
the generation function, which was parametrised as in equation (|8|), with pj r 9en = 0.8151 
and c 9en = 0. 

Using the most recent form factor parametrisation of equation (|Tl~|) the following results 
were obtained: 



A^l) | V cb | = (35.5 ± 1.4) x 1(T 3 
p\ = 1.34 ±0.14 
BR(B° d -> D*+l-Vf) = (4.70 ± 0.12)% 

where the last quantity was obtained by integrating the differential decay width. The 
correlation between the two fitted parameters was 0.94. Figure ^ shows the comparison 
between the real data and the result of the fit. 

It should be noted that the fit was performed separately on 1992-1993 and 1994-1995 
data sets, and then the results have been averaged. Individual results obtained with the 
two data sets are in agreement, as can be seen in Table |3|. 

The same table also contains the results obtained when using the other parametrisa- 
tions of form factors discussed previously. In detail, the Taylor expansion of equation 
(H) was employed, by assuming a linear expansion (c=0), by imposing the constraint of 
equation (0) for the curvature c or by fitting ftp and c as independent free parameters. In 
this last case, the correlation coefficients with |V^,|.F(1) were 0.82 and 0.71 respectively, 
the mutual correlation was 0.97. The last two entries in the table are reported to show the 



consistency with the published result of reference |18| . The measurement was performed 
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Fit Method (sample) 


Kb A(l)xlO a 


pX 


Pi 


c 


Eq. (11) 92-93 
Eq. (11) 94-95 
Eq. (|11|) 92-95 
Ref. @ 92-95 
Eq. @ 92-95 
Eq. © 92-95 
Linear 92-95 
Linear 94 only 
Ref. p| 


35.8 ±2.5 
35.2 ± 1.8 

o r r i i /i 

35.5 ± 1.4 

35.9 ± 1.6 

35.8 ± 1.4 

36.9 ± 1.9 

34.6 ± 1.3 
36.4 ± 1.5 
35.9 ±2.2 


1.30 ±0.24 
1.40 ±0.18 

"1 O A i r\ 1 A 

1.34 ± 0.14 
-0.0009 ±0.021 


1.22 ±0.14 
1.59 ±0.41 
0.90 ±0.10 
0.84 ±0.12 
0.74 ±0.20 


0.66p^-0.11 
1.4 ±0.9 






Table 3: Results of different fits to the experimental data. Results in the fourth line are 
obtained assuming the form factor representation of reference || , where the expansion is 
performed directly on z. For this reason the first order Taylor coefficient (p\ ) cannot be 
compared directly to the corresponding one of reference [JTTJ in the first three lines. 



by using the same data sample (1994) and the same model assumptions for the signal 
and the background as in that previous publication, but applying the new data selection 
and w reconstruction. 



6 Systematic Uncertainties 

The individual sources of systematic errors are reported in Table ^ and are described 
in detail below. Uncertainties in the overall normalisation, in the knowledge of the 
selection efficiency and of the composition of the sample, including the modelling of 
the background, and about the agreement between the experimental and the simulated 
resolution may affect the results. They were all considered as sources of systematic error. 

The fit was performed several times, by varying in turn all the parameters which deter- 
mine the normalisation (see equation (|8|)) within their allowed range. The corresponding 
variations of the measured quantities were added quadratically to the systematic error. 

The efficiency depends on the detector performance in track reconstruction, lepton 
identification and secondary vertex reconstruction. Tracks from charged particles may 
be lost because of cracks in the tracking device, or because of hard scattering of the 
particle by the detector frames. Electrons and low momentum ir* are more sensitive to 
this last effect. A conservative error of ±1% per track was assigned, based on studies 
of the detector material (performed using electrons from photon conversion) and of the 
TPC cracks. 

The actual efficiency for lepton identification was measured exploiting samples of real 
data tagged independently. Muons from r decays and from the process 77 — > /x + \i~ were 
used to explore all the relevant kinematic range. The values of the experimental and 
predicted efficiencies were consistent within ±2%. Electrons from photon conversion and 
from the radiative Bhabha process were also used. Compared to the simulation, a relative 
efficiency of (94 ±2)% was found, where the error is due to the systematic difference 
between the two samples. This ratio does not depend on the particle momentum. 

To provide an accurate description of the algorithm employed for vertex reconstruction, 
the simulation was tuned following the procedure of reference J34], developed for the 



precise measurement of Rb. The efficiency was then determined by comparing in the real 
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data and in the simulation the fraction of vertices reconstructed in a sample with high 
momentum leptons. The ratio between the efficiency found in the real data to the one 
found in the simulation was 1.01±0.01. The average number of charged particles forming 
the inclusive vertex in the simulation was slightly larger than in the real data. This was 
attributed to a small loss in the efficiency to assign the charged particle tracks to the 
secondary vertex. The ratio of the vertex reconstruction efficiencies was estimated to be 
0.99 ±0.01. 

Because of the cuts on the lepton momentum and decay length significance Sc (see 
Section |4.1|) , the efficiency depends on the average fraction of the beam energy ac- 
tually carried by the B® meson, < xe >■ Events were generated assuming the Pe- 
terson fragmentation function ||35| , tuned so as to reproduce the measured value of 
< xe > = 0.702 ± 0.008 |p6| ; they were then reweighted in the fit in order to allow 
for a variation of ± 0.008 in < xe >, and the consequent change of the fitted parameters 
was propagated into the errors. 

Model dependent uncertainties may be introduced by the kinematic cuts on p\ and 
fi 2 as well. They were determined following the iterative procedure applied in refer- 
ence |18| . The simulated spectrum was corrected to the measured values and the effi- 
ciency computed again. The efficiency varied by about ± 1%. This was taken as the 
systematic error, and no further iteration was performed. 

Each source of background was changed by its error as given in Table ^| and the 
variation of the results was propagated into the error. 

The reconstruction efficiency and the w resolution depend on the multiplicity of 
charged particles in D° decays, improving with higher values. It should be noted that 
zero-prong D° decays are also collected, albeit with smaller efficiency and with worse w 
resolution, in all the cases in which a secondary vertex can be formed by the lepton and 
the 7r* alone. The simulation was tuned to the results of the MARK III measurement 
|55fl . The relative fractions of events with 0, 2, 4, >4 charged prongs were then varied 
within their errors to compute the systematic error. 

The fraction of K° produced per D° decay was varied as well, to account for possible 
loss of efficiency and degradation of resolution due to the presence of K\. 

The systematic error due to the knowledge of the D** spectra was determined following 



the prescriptions of the LEP | V^,|Working Group, as explained in detail in section [572 



All quantities relevant to the determination of the w resolution were studied. Agree- 
ment was found between the distributions of the D* + energy in the real data and in the 
simulation; the angular resolution on the D* + direction was checked by comparing the 
relative angle between the n* and the D° directions and again very good agreement was 
found. 

The estimate of the B® energy depends on the hermeticity of the detector. To ver- 
ify that cracks were properly simulated, a sample of b enriched events was provided by 
6-tagging one hemisphere and analysing the other (unbiased). Only hemispheres with- 
out identified leptons were considered, in order to avoid possible distorsions due to the 
presence of a neutrino. The procedure was applied to the experimental data and to the 
simulation, and the visible energy was compared in the two cases. Depending on the year, 
the energy seen in the real data (about 37 GeV) exceeded the one predicted by about 
200-400 MeV. The main source of discrepancy was attributed to the tracking. Due to the 
smallness of the effect, no further investigation was performed but two different correction 
procedures were followed: either the visible energy in the simulation was increased by the 
relevant amount, or else a correction was computed depending on the fraction of charged 
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energy seen in the event. The systematic error was chosen as the maximum difference 
between the result so obtained and the one without fine tuning. 

The B® angular resolution was compared in the real data and in the simulation by 
inspecting the angle between the B® and the D* + i~ directions. The RMS widths of the 
two distributions were identical within errors (42.5±0.6 and 42. 6± 0.6 mrad respectively). 
The systematic error on \ V c b\ was computed by repeating the fit without the fi 2 constraint. 
The w resolution depends on the precision of the vertex reconstruction, which improves 
when more charged tracks form the vertex. Events in the simulation were rescaled to 
correct for the small discrepancy in the vertex multiplicity discussed previously and the 
fit was repeated. The difference was negligible. 

The systematic error induced by the fitting procedure was determined by varying the 
number and the size of the bins, and by removing the (few) events outside the physically 
allowed w region. The effect of the other cuts applied in the analysis was checked by 
varying them in the ranges : 

• the Am cut between 0.15 and 0.20 GeV/c 2 , 

• the p[ cut between 0.8 and 1.25 GeV/c, 

• the /i 2 cut between 0. and 5. GeV 2 /c 4 

• the Sc cut between 2.5 and 6.5 

The efficiency and purity of the sample vary by more than 50% in these ranges, and 
most of the induced variations are compatible with statistical fluctuations. They were 
conservatively assumed as systematic errors. As a further check, the analysis was per- 
formed separately for electrons and muons. Excellent agreement was found. All the errors 
were added in quadrature to determine the final systematic uncertainty. 

7 Extraction of the form factor 

The result of the previous section was obtained in the framework of a specific model. It 
is in principle possible to extract the differential decay width, dT/dw, from the experimen- 
tal data. To cope with the non- negligible smearing due to the experimental resolution, 
an unfolding procedure was applied |37| , p8| . With this same technique the Isgur Wise 
function, the universal form factor expected in the framework of the HEQT, was also 
extracted. 

The simulated events which survived the selection were first grouped in ten bins, ac- 
cording to the value of w gen at generation. Because of the finite experimental resolution, 
events lying inside a given bin in w gen populated several bins in the reconstructed w rec 
distribution. For each w gen bin, a corresponding w rec histogram was obtained. To over- 
constrain the fit, the new histograms consisted of twelve bins. The linear combination of 
these ten histograms was fitted to the real data distribution. The ten parameters of the 
fit determined the normalisation coefficients for each simulation histogram. The unfolded 
differential decay width was finally obtained by binning the simulated events according 
to the value of w gen and scaling the resulting histograms with the fitted parameters. 

To avoid spurious bin-to-bin oscillations, typical of such an unfolding method, a reg- 
ularisation term was added to the x 2 ; which is proportional to the second derivative of 
the unfolded results: 

X 2 reg = r • E^K/m - /*) - (fi ~ /*-i)| 2 T-J \f"(x)fdx (12) 

The regularisation parameter r is in principle arbitrary. Too small values lead to oscil- 
lating solutions, whereas large values produce flat solutions with small errors and strong 
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positive correlations among parameters. Several fits were performed with r ranging from 
0.01 to 1.0. To test the method, the unfolded distributions were fitted with the function 
of equation ( jTT| ) neglecting bin-to- bin correlations. The values obtained for Ay \V c b\ and 
p\ were always well compatible with those given in section [5], but their errors depend 
on the choice of r (lower values leading to higher errors). Choosing r=0.20, the same 
errors as the ones of section |5| were obtained. The corresponding unfolded spectra are 
presented in Figures |5|(a,b); the dashed curve overlayed shows the result obtained when 
fitting neglecting bin-to-bin correlation. To remove the sensitivity of the errors to the 
choice of t, fits were finally performed properly accounting for bin-to-bin correlations: 
they are represented by the continuous line in Figure |5](b). The result was: 

A x -\V ch \ = (36.1 ± 1.4) x 10" 3 
p\ = 1.38 ±0.15 

independent of the choice of r. The small difference from the values presented in section 
|B~3| is interpreted as the systematic error due to the unfolding procedure. The unfolded 
data and their error matrix are presented in Table [5[ 



8 Conclusions 

A sample of about 5000 B® — > D* + £~u e decays was obtained by means of the method 
of the inclusive 7r* tagging, originally developed at LEP by the DELPHI collaboration. 
The use of the large data set, and the excellent detector performance allowed the precise 
measurement of the product |Vd,|-.Ai(l) and of the B® "radius" p\ x , following the most 



recent parametrisation of the Isgur-Wise function proposed in reference 1 1 



| V ch | = (35.5 ± 1.4(stat.)+| 3 (syst.)) x 10" 3 

p\ = 1.34±0.14(stat.)+;|!(syst.) 
BR{B° d -> D* + £-v e ) = (4.70±0.13±°; 3 ?)% 

Using the value .Ai(l) ~ ^"(1) = 0.91 ±0.03, the following value of \V c b\ is obtained: 

| V cb | = (39.0 ± 1.5(stat.)+^(syst. exp.) ± 1.3(syst. th)) x 10~ 3 

These results agree with the present world average (see reference f25|| ). They supersede 
the previous DELPHI measurement of reference 
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Table 4: Contributions to the systematic uncertainties. The values used for the parame- 
ters relevant to this analysis are reported in the second column. Errors for \V c b\ and BR 
(B® — > D* + £~v e ) are relative and given in %; the errors for p^are absolute. 
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Table 5: The unfolded differential decay width. First line: w bin centre. Second line: 
partial decay width in the corresponding bin, divided by the total width, expressed in 
percent. Other lines: corresponding error matrix expressed in permill. All the values 
were obtained with the regularisation constant r = 0.20. 
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Figure 1: Mass difference Md-k — Md- 

a) opposite charge, real data (dots with error bars), same charge (shaded area), normalised 
to the side band defined in the text. The D* + signal is clearly visible. 

b) opposite charge, real data after subtraction of the combinatorial background (dots 
with error bars). This agrees well with the resonant contribution from simulation (shaded 
area) ; 

c) simulation: combinatorial background from opposite charge (dots), which is consistent 
with the same charge combinations normalised in the side band (shaded area). 
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Figure 2: w resolution. Upper plot : 1992-1993 analysis; lower plot: 1994-1995 analysis. 
Dots: experimental resolution without exploiting kinematic constraints. Since 1994 
three-dimensional vertex reconstruction helped improve the resolution. 
Shaded area: further improvement due to the requirement /i 2 = (see text). 
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Figure 3: Squared missing mass distributions. Real data (dots with error bars) after 
subtraction of all the background apart from D** are compared to the sum of the D* + 
and Z}**contributions as predicted by the simulation; the dark area represents the D**. 
The vertical line shows the position of the cut. 
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Figure 4: Fit to the w distributions. Dots with error bars: data; dark shaded area: com- 
binatorial background; light shaded area: other backgrounds, including D**; histogram: 
all components, the unshaded area corresponds to the decay B® — ► D* + £~u e . 
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Figure 5: Unfolded distributions in the real data. Upper plot: differential decay width. 
Lower plot: decay form factor as in equation fllCf). The dotted line shows the results of 
a fit to the histograms, neglecting bin-to-bin correlations. The continuous line shows the 
result obtained when including the statistical correlations among the bins. 



